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ABSTRACT 

We report the discovery of a 38.5 ms X-ray pulsar in observations of the soft 7-ray source 
IGR J18490-0000 with the Rossi X-ray Timing Explorer (RXTE). PSR J1849-0001 is spinning down 
rapidly with period derivative 1.42 x 10~ 14 s s" 1 , yielding a spin-down luminosity E = 9.8 x 10 36 
erg s -1 , characteristic age t c = P/2P = 42.9 kyr, and surface dipole magnetic field strength B s = 
7.5 x 10 11 G. Within the INTEGRAL/IBIS error circle lies a point-like XMM-Newton and Chandra X- 
ray source that shows evidence of faint extended emission consistent with a pulsar wind nebula (PWN) . 
The XMM-Newton spectrum of the point source is well fitted by an absorbed power-law model with 
photon index T PSR = 1.1±0.2, iV H = (4.3± 0.6) x 10 22 cm" 2 , and F PSR (2-10 keV) = (3.8±0.3) x KT 12 
erg cm -2 s -1 , while the spectral parameters of the extended emission are roughly Tpwn ~ 2.1 and 
F PWN (2 - lOkeV) « 9 x KT 13 erg cm" 2 s -1 . IGR J18490-0000 is also coincident with the compact 
TeV source HESS J1849— 000. For an assumed distance of 7 kpc in the Scutum arm tangent region, 
the 0.35 - 10 TeV luminosity of HESS J1849-000 is 0.13% of the pulsar's spin-down energy while the 
ratio -F(0.35 — 10 TcV)/Fpwn(2 — lOkeV) « 2. These properties are consistent with leptonic models 
of TeV emission from PWNe, with PSR J1849— 0001 in a stage of transition from a synchrotron X-ray 
source to an inverse Compton 7-ray source. 

Subject headings: ISM: supernova remnants pulsars: individual (HESS J1849— 000, 

IGR J18490-0000, PSR J1849-0001) — stars: neutron 



1. INTRODUCTION 

The detection of 10 12 eV emission associated with su- 
pernova products in the Galaxy has opened up a new 
window on the evolution of these energetic stellar rem- 
nants. More than 2/3 of the now > 60 Galactic TeV 
sources are supernova remnants or pulsar wind nebu- 
lae (PWNe), the latter being the largest clasfl Many 
of these TeV PWNe are spatially offset from middle- 
aged (t c ~ 10 4 — 10 5 year old) pulsars, and are of- 
ten more extended and more luminous than their X-ray 
PWNe counterparts. This TeV emission is likely to be 
inverse Compton scattered radiation from relic electrons 
produced b y the pulsar in an earlier stage of energetic 
spin-down dde Jager fc Diannati-Atal f2008; Zh ang et al.l 
120081: iMattana et al.ll2009al) In contrast, younger (r c ~ 
10 3 year old) pulsars are associated with compact TeV 
sources that are co-located with their X-ray PWNe. In 
these younger systems the high magnetic fields make for 
efficient synchrotron X-ray sources, and inefficient in- 
verse Compton TeV emission. 

HESS J1849— 000 is one of the fainter sources detected 
in the HESS Galactic Plane survey (EE aronian et al.1 
120051 l2006f ) and subsequent dedicated observations, with 
significance of 6.4cr, flux > 350 G eV of ~ 15 mCrab , 
and only weak indication of extent (jTerrier et al.| [2008). 
HESS J1849— 000 is coincident in position with the 
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INTEGRAL/IBIS source IGR J18490-0000 that was 
discovered during a survey of the tange nt regions of 
the Sagittarius and S cutum spiral arms (jMolkov et al.1 
I2004t iBird et aUfe OOr]!) . Follow-up X-ra y observation of 
IGR J18490-0000 by iRodriguez eFaLl ((20081) with the 
Swift X-ray telescope (XRT) located a highly absorbed 
X-ray point source within the INTEGRAL/IBIS error 
circle. These authors suggested a Galactic X-ray binary 
origin for IGR J18490-0000 based on a K-band Two 
Micro All Sky Survey star located within the Swift XRT 
error circle. However, that association was excluded by 
the precise X-ray position obtained in a br ief Chandra 
High R esolution Cam e ra (H RC) observation ([Ratti et al.l 
Mm . iTerrier et all l|2008l ) used imaging and spectro- 
scopic evidence from an XMM-Newton observation to 
argue that IGR J18490-0000 is instead a pulsar/PWN. 
As such, it would join the dozen hard X- ray members of 
this class detected by INT EGRAL/IBIS (|Mattana et all 
I2009bb iRenaud et al.ll2010ft . 

In Section 2 we present in detail the XMM-Newton 
imaging and spectral data that support a pulsar/PWN 
interpretation for IGR J18490-0000. In Section 3, we re- 
port the discovery using Rossi X-ray Timming Explorer 
{RXTE) of PSR J1849-0001 and its spin-down, which 
verifies this conjecture. In Section 4, we discuss the prop- 
erties of HESS J1849-000 in the context of the spin- 
down parameters of PSR J1849-0001. 

2. XMM-Newton OBSERVATION 

An 11 ks XMM-Newton observation of 
IGR J18490-0000 (ObsID 0306170201) was ac- 
quired on 2006 April 3 using the Europea n Photon 
Imaging Camera (EPIC; iTurner et ail 120031 ). EPIC 
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Figure 1. XMM-Newton EPIC X-ray observation of the INTEGRAL/IBIS source ICR J18490-0000. Left: exposure corrected 
and smoothed EPIC MOS image scaled in intensity to highlight the faint diffuse emission near the point s ource. This source , 
the putative pulsar, lies within the INTEGRAL/IBIS r = I'A positional uncertainty of ICR J18490-0000 ijTerrier et alJl2008t 
dashed circle), as does the centroid of HESS J1849— 000, indicated by the triangle. Right: XMM-Newton EPIC pn (black) and 
EPIC MOS (red) spectra of the point source fitted to the model described in the text. The lower panel shows residuals from 
the best fit in units of la. 

consists of three sensors operating in parallel, one pn 
and two MOS CCD cameras. These instruments are 
sensitive to X-rays in the 0.2-12 keV range with energy 
resolution AE/E 



0.1/^/25(keV). All three cameras 
were operated in full-frame mode, using the thin and 
medium filters for the EPIC pn and MOS, respectively. 
The target was placed at the default EPIC pn focal 
plane location for a point source. The time resolution 
of 73.4 ms and 2.7 s for the EPIC pn and MOS, 
respectively, are insufficient to search for a typical pulsar 
signal. 

The data were processed using the SAS version xmm- 
sas_20061026_1802-6.6.0 pipeline, and were analyzed us- 
ing both the SAS and FTOOLS software packages. The 
observation was free of significant particle background 
contamination and provided a near continuous 11 ks of 
good observing time for the EPIC MOS and 9.9 ks for 
EPIC pn. 

Figure Q] (left) displays the 2-10 keV XMM-Newton 
EPIC MOS X-ray image centered on ICR J18490-0000. 
The image has been exposure corrected and smoothed 
using a a = 3'.'8 Gaussian kernel and scaled to highlight 
the faint diffuse emission near the central bright source 
at coordinates (J2000.0) R.A. = 18 h 49 m 01?62, decl. = 
— 00°01'17."7. This source, at the position of the previ- 
ously detected Swift XRT source, lies at the center of the 
1.'4 radius error circle of IGR J18490-0000. There are 
no other prominent X-ray sources in the full 30' diame- 
ter field of this instrument. Together with the positional 
coincidence, its flux and spectrum (see below) leave no 
doubt that XMMU J184901. 6-000117 is the counterpart 
of IGR J18490— 0000. Faint, diffuse emission surrounds 
the point source, and a prominent feature extends up 
to 1 to the southwest; these properties are suggestive 
of a PWN. A 1.2 ks Chandra HRC image obtained on 



2008 February 16 (|Ratti et alJl2010j ) shows that most of 
the XMM-Newton source flux is contained in a point-like 
component. 

2.1. XMM-Newton Point-source Spectrum 

A spectrum of XMMU J184901.6-000117 was ex- 
tracted from each EPIC camera using a 30" radius aper- 
ture centered on the point source; data from the two 
MOS cameras were combined into a single spectrum. The 
source background was estimated using counts from an 
60" aperture placed just south of the source region, on 
the same CCD of each camera. These source spectra were 
grouped with a minimum of 50 counts per spect ral chan- 
nel an d fitted using the XSPEC fitting package ([Arnaud I 
1996). EPIC pn and MOS spectra were fitted simulta- 
neously with a single absorbed power-law model, allow- 
ing independent flux normalization (Figure [TJ right); we 
report the average of the fluxes measured by the two 
EPIC cameras. In the following, all derived luminosities 
are corrected for interstellar absorption, fluxes are uncor- 
rected, and errors are at the 90% confidence level. The 
data yield a good fit with a reduced xl — 0-90 for 72 de- 
grees of freedom (dof). The hard spectrum has photon 



(4.3 ± 0.6) x 10 22 cnr 



(3.8±0.3)xl0" 12 erg cm" 2 s" 



index r PSR = 1.1 ±0.2, 
andF PSR (2-10keV) 

The Galactic coordinates of IGR J18490-0000, (£, b) 
(32.°64, +0.°53) are coincident with the tangent point of 
the Scutum arm. Considering this and the large fit- 
ted N K , we hypothesize that IGR J18490-0000 is lo- 
cated at a distance of 7 kpc appropriate for the Scu- 
tum tangent region rather than residing in the nearer 
Sagittarius arm. Its 2 — 10 keV X-ray luminosity is then 
L PS r(2 - lOkeV) = 2.8 x 10 34 d| erg s" 1 . 

We also fitted the XMM-Newton and INTE- 
GRAL /IBIS spectra simultaneously, finding that a sim- 
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Figure 2. Radial profile of the emission from IGR J18490-0000 
in the EPIC pn. The dotted line is the average King profile for a 
r = 1.1 power-law spectrum; the red points are background taken 
at a symmetric position with respect to the instrument optical 
axis. An additional extended component visible up to 150" from 
the point source is well fitted by a Gaussian profile of <r = 75" 
(dashed line). The solid line is the total fitted profile. 

pie power-law model gives a reasonable \ 2 but a rather 
poor fit to the INTEGRAL/IBIS data (observation and 
analysis of the IGR J18490— 0000 data as described in 
Terrier et al. 2008). An excellent result is obtained 
by fitting a broken power law (xl = 0.95 for 140 dof) 
with a break energy of 23 keV. The best fit indices of 
T = 1.2 ± 0.1 and T = 1.7 ± 0.2, for the low and high 
energies, respectively, are co nsistent with those obtained 
for the spectra individually (jTerrier et al.ll2008l ). 

2.2. XMM-Newton Extended Emission 

To quantify the spatial structure of the XMM-Newton 
source, we made a radial profile of the counts distribution 
centered on the point source and compared it with the 
point spread function (PSF). The latter was constructed 
using instrumental profiles at several energies weighted to 
match the measured spectral index. Figure [2] shows the 
measured EPIC pn profile compared to the average PSF 
for this spectrum and location in the field. Background 
has been taken from a symmetric position with respect 
to the optical axis. There is a clear excess of emission 
from 20" to 150" radius. We fitted the radial counts 
distribution to the sum of a King profile to represent the 
point source, and a Gaussian to represent the extended 
emission. A Gaussian component 150" ± 15" in diameter 
improves the quality of the fit by « 9c, proving that the 
source is surrounded by a nebula or a halo. 

We extracted the spectrum of the extended nebula in 
an annulus surrounding XMMU J184901.6-000117 from 
30" to 150". To account for contamination from the 
wings of the bright point source, we fitted the point 
source and extended emission spectra simultaneously. 
We used separate power laws but the same Ah for both 
spectra. We added a constant fraction of the point- 
source spectrum to the diffuse emission. This fraction 



was estimated from the PSF as 12% and 14% for EPIC 
pn and MOS, respectively, at a radius of 30". To ac- 
count for the uncertainty in the estimate of these frac- 
tions, we varied them from 7% to 19% of the point-source 
flux, adding the corresponding systematic uncertainty to 
errors on the fitted parameters. The resulting diffuse 
emission spectrum has Tpwn = 2.1 ± 0.3 and a flux of 
F PW n(2- lOkeV) = (9 ±2) x 10~ 13 erg cm" 2 s -1 . The 
reduced chi-squared is \1 — 0-90 for 311 dof. 

The extended nebula with a steeper spectrum suggests 
that IGR J18490-0000 is a pulsar/PWN system. Alter- 
natively, given the large Ah, some of the diffuse emis- 
sion co uld be a dust-scattered h alo of the bright point 
source (jPredehl fc Schmittl Il995[ k however, that could 
not account for its asymmetric extension. We note that 
the nebula is rather faint, m 25% of the point source 
flux, but it is still in the range observed for PWNe 
(jKargaltsev fc P avlovl [20081 ) . 

3. RXTE OBSERVATIONS 

To search for the expected pulsar signal from 
IGR J18490-0000, we obtained 112.3 ks of RXTE ob- 
servations, pointed at the Chandra source position, span- 
ning 2010 November 25 - December 15. The schedule 
was designed to obtain a phase-connected timing solu- 
tion that would measure both P and P. The data used 
here w ere collected with th e Proportional Counter Array 
(PCA; [Tahoda et al.lll996f ) in the GoodXenon mode with 
an average of 1.75 out of the five proportional counter 
units (PCUs) active. In this mode, photons are time- 
tagged to 0.9 /xs preci sion and have a n absolute time 
accuracy of < 100 /zs (jRots et all 119981 ). The effective 
area of five combined detectors is ~ 6500 cm 2 at 10 keV 
with a roughly circular field of view of ~ 1° FWHM. 
Spectral information is available in the 2—60 keV energy 
band with a resolution of ~ 16% at 6 keV. 

Standard time filtering was applied to the PCA 
production data, rejecting intervals of South Atlantic 
Anomaly passage, Earth occultation, and other periods 
of high particle activity. The photon arrival times were 
transformed to the solar system barycenter in Barycen- 
tric Dynamical Time (TDB) using the JPL DE200 
ephemeris and the coordinates given in Table [I] obtained 
from Chandra HRC observation ObsID 7398. These coor- 
dinates are derived from a centroid calculation of 17 pho- 
tons within a 2" radius aperture and have nominal un- 
certainty of 0"6. We note that these coordinates differ 
by 0"35 from those reported in Ratti et al. (2010). 

3.1. RXTE Timing Analysis 

We restricted the analysis to the 2 — 20 keV energy 
range (PCA channels 2 — 50) from the top Xenon layer of 
each PCU to optimize the signal-to-noise. A fast Fourier 
transform of the initial data set (ObsID 95309-0 TO T 
02; 7 ks) revealed a highly significant signal of period 
P = 38.52 ms. Subsequently, for each satellite orbit 
we extracted a pulse profile correspondin g to the pea k 
power as determined by the Rayleigh test (|Struttl ll880). 
also known as Z\ (jBuccheri et al.l fl983f ). The resulting 
profiles were cross-correlated, shifted, and summed to 
generate a master pulse profile template. Each individ- 
ual profile was then cross correlated with this template 
to determine its time of arrival (TOA) and uncertainty. 
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Figure 3. Discovery timing results for the X-ray pulsar PSR J1849— 0001 using RXTE data. Pulse phase residuals (left) and folded light 
curve (right) of PSR J1849— 0001 in the 2—20 keV band were generated using the quadratic ephemeris of Table [T] Phase zero of the light 
curve corresponds to the ephemeris epoch; two cycles are shown for clarity. 

tions. Similarly, standard PCA responses for each PCU 
were generated at each ObsID and averaged. In fitting 
the pulsed flux, the unpulsed emission provides a near 
perfect background estimate and was taken from the 0.4 
phase bins range corresponding to the region of flat min- 
imum in the pulse profile around phase zero. 

The merged RXTE spectrum was fitted in the 2 — 
20 keV range using a simple absorbed power-law model 
with the interstellar absorption held fixed at N-r — 
4.3 x 10 22 cm~ 2 determined from the XMM-Newton 
fit (see Figure 4). The resulting best- fit photon in- 



Parameter Value 

R.A. (J2000.0) a 18 h 49 m 01. s 61 

Decl. (J2000.0) a -00°01'17."6 

Epoch (MJD TDB) 55535.285052933 

Period b , P 38.518943432(12) ms 

Period derivative 11 , P 1.4224(58) X 10 -14 s s" 1 

Range of dates (MJD) 55525.6 - 55545.2 

Spin-down luminosity, E 9.8 X 10 36 erg s — 1 

Characteristic age, t c 42.9 kyr 

Surface dipole magnetic field, B s 7.5 X 10 11 G 

a Chandra position with a nominal uncertainty of 0'.'6 (see the 
text). 

b TEMPO lcr uncertainties given in parentheses. 

These TOAs were then iteratively fitted to a quadratic 
ephemeris. 

The resulting unique ephemeris is presented in Ta- 
ble [1] and the phase residuals are shown in Figure El 
The residuals are all less than 0.05 cycles and do not 
require a higher derivative. The derived properties 
of PSR J1849— 0001 are: spin-down luminosity E — 
4n 2 IP/P 3 — 9.8 x 10 36 erg s _1 , characteristic age r c = 
P/2P = 42.9 kyr, and surface dipole magnetic field 
strength B s = 3.2 x 10 19 (PP) 1 / 2 = 7.5 x 10 11 G. 

Figure [3] displays the summed pulse profile using all 
the 2—20 keV data folded on the final ephemeris. It has 
a symmetric, single-peaked structure. The raw pulsed 
fraction is 3.8%, uncorrected for PCA instrument or as- 
trophysical background, difficult to determine with suffi- 
cient accuracy for this purposes. We see no energy depen- 
dence of the pulse profile when subdividing the 2-20 keV 
band. 



3.2. RXTE Spectral Analysis 
The spectrum of the pulsed flux from PSR J1849- 



0001 



can be isolated by phase-resolved spectroscopy. We used 
the fasebin software to construct phase-dependent spec- 
tra based on the ephemeris of Table [T] For each ObsID 
and PCU we constructed spectra from counts detected in 
the top Xenon layer only and combined them to produce 
a single spectrum per PCU for the entire set of observa- 



dex is T = 1.31^3. The 2 - 10 kcV pulsed flux is 
8.9 x 10~ 13 erg cm -2 s _1 , which represents ~ 25% of 
the point-source flux of PSR J1849— 0001 measured by 
XMM-Newton. The corresponding pulsed luminosity (as- 
sumed isotropic) is 6.7 x 10 33 dj erg s _1 . The reduced 
chi-squared is xl = 0.70 for 29 dof. 

As expected, the RXTE measured spectral slope 
is intermediate between the XMM-Newton and INTE- 
GRAL/IBIS values, indicating a steepening around 10 — 
20 keV. 

4. DISCUSSION AND CONCLUSIONS 

The discovery of PSR J1849-0001 u sing RXTE 
verifies the conjecture of iTerrier et al.l ((2008) that 
XMMU J184901.6-000117, IGR J18490-0000, and 
HESS J1849— 000 are all manifestations of a young pul- 
sar /PWN system. Even though the RXTE field of 
view is wider than those of XMM-Newton and Chan- 
dra, there is little doubt that PSR J1849-0001 is 
the compact source in IGR J18490— 0000, given the 
morphological and spectral evidence, and the com- 
patibility of the RXTE measured puls ed flux with 
that of XMMU J184901.6-000117. ITerrier et all 
(2008) estimated that the spin-down luminosity of 
PSR J1849-0001 would be 9 x 10 36 erg s" 1 based on 
the em pirical correlation between E and Tpsr of lGottheTfl 
(2003); this turns out to have been an accurate predic- 
tion. 

Assuming a distance of 7 kpc as proposed in Sec- 
tion 2.1, the X-ray luminosities of the pulsar and 
PWN, Lp S r(2 - lOkeV) = 2.9 x lO^d 2 . and 
Lpwn(2 — lOkeV) = 6.6 x \Q~ A Ed 2 -, respectively, 
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Figure 4. RXTE spectrum of pulsed flux from PSR J1849-0001 
obtained by subtracting the off-peak spectrum from the on-peak 
spectrum, and fitting to an absorbed power-law model (see the 
text). 



are consistent with the ran ge of pulsar/PWN systems 
(|Kargaltsev fc Pavlovl [200l. The 20 - 100 keV flux 
of IGR J18490-0000 measured by INTE GRAL /IBIS 
(2 x 10~ n erg cm~ 2 s -1 : lTerrier e t al. 200§) corresponds 
to L(20 - 100 keV) = 0.012 Ed%. Given the 2-20 keV 
flux and spectral information from XMM-Newton and 
RXTE, the INTEGRAL/IBIS source is likely to be dom- 
inated by the pulsar rather than the PWN. 

The TeV flux from HESS J1849- 000 (2.2 x 10~ 12 
erg cm~ 2 s _1 ; iTerrier et al.l I2008T) corresponds to 

L(0.35 - 10 TeV) = 1.3 x lO" 3 ^?.. Such a small ef- 
ficiency of converting spin-down power to high-energy 
radiation is typical of high E pulsars. The ratio F(0.35 — 
10TeV)/F PWN (2 - lOkcV) 2 falls squarely on the in- 
verse cor relation between this q uantity and E that was 
fitted bv iMattana et all (|200 9a). and modeled by them 
in terms of an evolving PWN emitting synchrotron X- 
rays and inverse Compton 7-rays. PSR J1849— 0001 is in 
transition from a synchrotron dominated X-ray PWN to 



an inverse Compton scattered TeV nebula, an expected 
phase through which a middle-aged pulsar will pass. 

We thank the RXTE project for making the time avail- 
able for this program, and the mission planners for care- 
fully scheduling the observations. This investigation is 
also based on observations obtained with XMM-Newton, 
an ESA science mission with instruments and contribu- 
tions directly funded by ESA Member States, and NASA. 
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